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Q ■ Abstract 

_J I Derivatives of fractional order with respect to time describe long-term memory effects. 

Using nonlinear differential equation with Caputo fractional derivative of arbitrary order 
a > 0, we obtain discrete maps with power-law memory. These maps are generalizations 
►-j^ , of well-known universal map. The memory in these maps means that their present state is 

_J ■ determined by all past states with power-law forms of weights. Discrete map equations are 

^T \ obtained by using the equivalence of the Cauchy-type problem for fractional differential 

f~^ ■ equation and the nonlinear Volterra integral equation of the second kind. 

o 

1 Introduction 

c^ ; Discrete maps are used for the study of evolution problems, possibly as a substitute of differ- 

ential equations [H El [3] . They lead to a much simpler formalism, which is particularly useful 
in numerical simulations. The universal discrete map is one of the most widely studied maps. 
It is a very important step in understanding the qualitative behavior of a wide class of sys- 
tems described by differential equations. The derivatives of noninteger orders [U Ej E] are a 
natural generalization of the ordinary differentiation of integer order. Fractional differentiation 
with respect to time is characterized by power-law memory effects. The discrete maps with 
memory are considered in Refs. [TJ IH1 El [ini [HI [El [13]. It is important to connect fractional 



differential equations and discrete maps with memory. In Ref. [13], we prove that the discrete 
maps with memory can be derived from differential equations with fractional derivatives. The 
fractional generalization of the universal map was derived [13] from a differential equation with 
Riemann-Liouville fractional derivatives. The Riemann-Liouville derivative has some notable 
disadvantages in physical applications such as the hyper- singular improper integral, where the 
order of singularity is higher than the dimension, and nonzero of the fractional derivative of 
constants, which would entail that dissipation does not vanish for a system in equilibrium. The 
desire to formulate initial value problems for physical systems leads to the use of Caputo frac- 
tional derivatives rather than Riemann-Liouville fractional derivatives. In this paper, we obtain 
a discrete map with memory from differential equations with Caputo fractional derivative of 
arbitrary order a > 0. The universal map with power-law memory is obtained by using the 
equivalence of the fractional differential equation and the Volterra integral equation. We reduce 
the Cauchy-type problem for the differential equations with the Caputo fractional derivative to 
a nonlinear Volterra integral equation of the second kind. The equivalence of this Cauchy-type 
problem and the correspondent Volterra equation was proved by Kilbas and Marzan in [15] [T6] . 
In Section 2, differential equations with integer derivative and universal maps without mem- 
ory are considered to fix notations and provide convenient references. In Section 3, fractional 
differential equations with Caputo derivative and correspondent discrete maps with memory 
are considered. A fractional generalization of the universal map is obtained from kicked differ- 
ential equations with the Caputo fractional derivative of arbitrary order a > 0. Finally, a short 
conclusion is given in section 4. 

2 Universal map without memory 

In this section, differential equations with derivative of second order and the universal map 
without memory are considered to fix notations and provide convenient references. 
Let us consider the equation of motion. 



oo 



in which perturbation is a periodic sequence of delta-function-type pulses (kicks) following with 
period T = 27r/z/, K is an amplitude of the pulses, D^ = d'^/dt'^, and G[x] is some real-valued 
function. It is well known that this differential equation can be represented in the form of the 
discrete map, 

Xn+l - Xn = Pn+lT, Pn+1 " Pn = -KT G[Xn]. (2) 

Equations ([2]) are called the universal map. For details, see for example [Il[2l[3]. 

Traditional method of derivation of the universal map equations from the differential equa- 
tions is considered in Section 5.1 of Ref. 0. We use another method of derivation of these 
equations to fix notations and provide convenient references. It is easy to obtain the universal 
map by using the equivalence of the differential equation and the Volterra integral equation. 
The Cauchy-type problem for the differential equations, 

Dixit) =p{t), (3) 

oo 

Dlpit) = -KG[x{t)]J2s{^ - k), (4) 

fc=i 

with the initial conditions 

x(0) = xo, p{0)=po (5) 

is equivalent to the universal map equations of the form 

n 

Xn+l = Xo+ Po{n + 1)T - KT^ ^ G[xk] (n+l-k), (6) 



n 



Pn+1 = PO 



KTj2G[xk]. (7) 



k=l 



To prove this statement we consider the nonlinear differential equation ([T]) on a finite interval 
[0, tf] of the real axis, with the initial conditions ([5]). The Cauchy-type problem of the form ([1]) 
and ([SD is equivalent to the Volterra integral equation, 

oo „i 

xit)=xo+Pot-KY, dTG[xiT)]6(^^-k^{t-T). (8) 

k=l ^^ 



For riT < t < {n + 1)T, we obtain 

n 

x{t) =Xo+ Pot -KTJ2 G[x{kT)] {t - kT). (9) 

fc=i 

Equations (jH]) and (^ give 

n 

p{t)=po-KTj2G[x{kT)]. (10) 



k=l 



The solution of the left side of the {n + l)th kick 

Xn+i = x{tn+i - 0) = lim x{T{n + 1) - e), (11) 

p„+i = p{tn+i - 0) = lim p(T(n + 1) - e), (12) 

e-*-0+ 

where t„_|_i = (n + 1)T gives the map equations (jS]) and (jTj). This ends the proof. Note that 
equations ([6]) and (jTj) can be rewritten in the form ([2]). Using equations iQ and ([7]), the 
differences x^+i — Xn and pri+i — P„ give equations ([2]) of the universal map. 
We note that equations ([2]) with G[x] = —x give the Anosov-type system 

Xn+l ~ Xn = Pn+lT, Pn+1 ~ Pn = KTXn- (13) 

If G[x] = sin(x), then equations (j2]) are 

Xn+l -Xn= Pn+lT, Pn+l - Pn = -KT sin(x„). (14) 

This map is known as the standard or Chirikov map [1]. 

3 Fractional equation and universal map with memory 



In Ref. [13] we consider nonlinear differential equations with Riemann-Liouville fractional 
derivatives. The discrete maps with memory are obtained from these equations. The Riemann- 
Liouville fractional derivative has some notable disadvantages in physical applications such as 
the hypersingular improper integral, where the order of singularity is higher than the dimension, 
and nonzero of the fractional derivative of constants, which would entail that dissipation does 



not vanish for a system in equilibrium. The desire to formulate initial value problems for 
physical systems leads to the use of Caputo fractional derivatives |H [5] rather than Riemann- 
Liouville fractional derivative. 

The left-sided Caputo fractional derivative [T71 UHl [191 S] of order a > is defined by 

where m — 1 < a < m and qI^ is the left-sided Riemann-Liouville fractional integral of order 
a > 0, that is defined by 

o^r/w = f^ r 77^^^' (^>o)- (16) 

r(") Jo [t-ry " 

The Caputo fractional derivative first computes an ordinary derivative followed by a fractional 
integral to achieve the desire order of fractional derivative. The Riemann-Liouville fractional 
derivative o-Df is computed in the reverse order. Integration by part of (IT5l) gives 

m-l i^_a 

The second term in Eq. (1171) regularizes the Caputo fractional derivative to avoid the potentially 
divergence from singular integration at t = 0. In addition, the Caputo fractional differentiation 
of a constant results in zero ^DfC = 0. The Riemann-Liouville fractional derivative of a 
constant need not be zero P]. 

If the Caputo fractional derivative is used instead of the Riemann-Liouville fractional deriva- 
tive, then the initial conditions for fractional dynamical systems are the same as those for the 
usual dynamical systems. The Caputo fractional derivatives can be more applicable to dynami- 
cal systems than the Riemann-Liouville derivatives. Note that the Caputo fractional derivatives 
can be used to formulate a self-consisted fractional vector calculus [T4] . 

We consider the nonlinear differential equation of order a, where 0<m — l<a<m, 

^D^x{t) = G[t,x{t)], (0<t<t/), (18) 

involving the Caputo fractional derivative ^-Df on a finite interval [0,t/] of the real axis, with 
the initial conditions 

{D^x){0) = Ck, k = 0,...,m-l. (19) 



Kilbas and Marzan [151 HB] proved the equivalence of the Cauchy-type problem of the form 
flTSj) . flTOl) and the Volterra integral equation of second kind 

m— 1 -. „( 

^ W = E 1^' + rv^ / ^^ ^[^' ^(^)] (^ - ^)""' (20) 

fc=o ■ \ ) Jo 

in the space C"""i[0,t/]. 

The basic theorem regarding the nonlinear differential equation involving the Caputo frac- 
tional derivative states that the Cauchy-type problem (ITSll . (fT9|) and the nonlinear Volterra 
integral equation ( l20l) are equivalent in the sense that, if x{t) G C[0,t/] satisfies one of these 
relations, then it also satisfies the other. In [151 IB] (see also [1], Theorem 3.24.) this theorem 
is proved by assuming that a function G[t,x] for any x G W C M. belong to C^(0,t/) with 
< 7 < 1, 7 < a. Here (7-^(0, t/) is the weighted space of functions f[t] given on (0,t/], such 
that r/[t] GC(0,t/). 

Let us consider a generalization of equation ([T]) in the form of the fractional differential 
equation, 

oo 

^D^x{t)+KG[x{t)]^5(--k) =0, (m-l<a<m), (21) 

fc=i 

where ^Df is the Caputo fractional derivative, with the initial conditions 

D'tx{0) = x'-q'^ (s = 0,l,...,m-l). (22) 

Using x^^\t) = Dfx{t), s = 0, 1, ...,m — 1, equation fl2T|) can be rewritten in the Hamilton form. 

Theorem. The Cauchy-type problem for the fractional differential equations 

Dlx^'\t)=x^'+^\t), (s = 0,l,...,m-2) (23) 

oo 

Cjj»-m+i^(m-i)^^^ = -KG[x{t)]J2s{:^ -k), (m - 1< a < m), (24) 



with the initial conditions 

x^'\0) = x'q'\ (s = 0, l,...,m-l) (25) 



is equivalent to the discrete map equations, 

'1 (k+s) 






1 = E ^(^ + 1)'^' - fT^^ E (^ + 1 - ^)°'''^G'N, (26) 



A;=0 ^ ' k=l 

Proof. Using the Kilbas-Marzan result for equation flTSjl with the function 

fe=i 

we obtain that the Cauchy-type problem ( 12T1) and ( l22l) is equivalent to the Volterra integral 



oo 

G[t,xit)] = -KG[xit)]Y,s{j; - k); 



(27) 



equation of second kind, 

m-l (fc) ^ oo „t 

-W = E^^'-f7^E/ ^r(t-r)"-G[x(r)]5(^-^), 

A:=0 ■ ^ '^ fc=l "^° 

in the space of continuously differentiable functions x{t) G C"^^^[0,tf]. 
If nT < t < (n + 1)T, then equation (l27Il gives 

m-l (fc) ^.^ m 

^(^) = E ir^' " ff^ E (^ - ^^)""' G[x(A:T)]. (28) 

fc=0 ■ ^ ' k=l 

Using the variables (l23!l . equation ( 128|) gives 

^^^^('^)= E ^^'-rY^^7yE(^-^^)""'"^^w^^)]' (29) 

fc=0 ■ ^ ^ k=l 

where s = 0, 1, ...,m — 1, nT < t < (n + l)T, m — l<a<m and we use r(z) = {z — l)T{z~ 1). 
The solution of the left side of the {n + l)th kick (fTTl) and ( fT2l) can be represented by equations 
fl26|l . where we use the condition of continuity x'^(t„ + 0) = x'^(t„ — 0), s = 0, 1, ..., m — 2. 
This ends the proof. D 



Equations fl26|) define a generalization of the universal map. This map is derived from a 
fractional differential equation with Caputo derivatives without any approximations. The main 
property of the suggested map is a long-term memory that means that their present state 
depends on all past states with a power-law form of weights. 



If G[x\ = sin(a;), then equations (!26l) define a generalization of standard map. For G[x] = 
—X, we have Anosov-type system with memory. 

In the case ofl<a<2,?7i = 2, we have the following universal map with memory: 

Xn+i = xo+ Po{n + 1)T - — - V (n + 1 - k^-^Glxk], (30) 

Pn+1 =Po- 7^7 -T V (n + 1 - kr-'Gix,]. (31) 

where a:„ = Xn and Pn = Xn . li a = m = 2, then equations ( l26l) give the universal map of 
the form ([6]) and ([7]) that is equivalent to equations ([2]). As a result, the usual universal map 
is a special case of this universal map with memory. 

4 Conclusion 

Equations for discrete maps with memory are suggested. The maps with power-law memory 
describe fractional dynamics of complex physical systems. The suggested map with memory is 
a generalization of well-known universal map. These maps are equivalent to the correspondent 
fractional kicked differential equations. To derive the map equations an approximation for 
fractional derivatives is not used. We obtain a discrete map with memory from fractional 
differential equation by using the equivalence of the Cauchy-type problem and the nonlinear 
Volterra integral equation of the second kind. 

Fractional differentiation with respect to time is characterized by power-law memory effects 
that correspond to intrinsic dissipative processes in the physical systems. Therefore, the uni- 
versal maps with memory have regular and strange attractors for some values of parameters K 
and a. The suggested universal maps with memory demonstrate a chaotic behavior with a new 
type of attractors. Numerical simulations of the universal map with memory prove that the 
nonlinear dynamical systems, which are described by the equations with fractional derivatives, 
exhibit a new type of chaotic motion. For some regions of parameters K and a these universal 
maps with memory demonstrate a new type of regular and strange attractors. The univer- 
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sal maps with power-law memory can be used to describe properties of regular and strange 
attractors of the fractional differential equations with kicks. 
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